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Abstract 

We study doubly polarised asymmetries of cc and bb mesons in 
hadro- and photo-production at low transverse momentum in non- 
relativistic QCD to high orders in the relative velocity of the pair, v. 
We give the complete set of expressions required for the asymmetries 
up to order a®. The asymmetries in the production of states are 
a stable measure of the polarised gluon densities. The asymmetries 
Xc,fe) and the various T states are stringent tests of the 

NRQCD scaling relations. 
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1 Introduction 


Recent progress in the understanding of cross sections for production of 
heavy quarkonium resonances has come through the non-relativistic QCD 
(NRQCD) reformulation of this problem |^. Factorisation of the physics at 
a short distance scale of order 1/m (where m is the mass of heavy quark) 
from long distance scales of order 1/mv and 1/mn^ (where v is the velocity 
of either of the heavy quark in the rest frame of the pair) has been proven 
in the NRQCD formalism for processes dominated by a large transverse mo¬ 
mentum P]. The resulting cross sections are a double power series in and 
the strong coupling as{m?)- For charmonium states, a numerical coincidence, 
~ as(m^), complicates the double expansion. 

The formalism has been successfully applied to large transverse momen¬ 
tum processes [|]. Interestingly, inclusive production cross sections for char¬ 
monium at low energies, dominated by low transverse momenta, also seem to 
have a good phenomenological description in terms of this approach iSB- 
The remaining phenomenological problem is to explain the X 1 /X 2 ratio ob¬ 
served in hadro-production. xq ,2 are hrst produced at order whereas Xi 
begins only at order OgU®. It was argued in p, |[ that a better understanding 
of this ratio may require higher order terms in ag and In the NRQCD 
series was evaluated to order In the dominant gluon fusion channel a 

previously missed term at order was found for xo ,2 and up to eleven 
terms at order were obtained for all y states. 

In this paper we extend the computations of 0] to polarisation asymme¬ 
tries, 

^ ^ ^ Ew hh'a{h, h') 

^ Ew^x(h,hO ’ ^ ^ 

where h and h' are the helicities of the beam and target respectively, and 
(j{h, h') is the cross section for hxed initial helicities. This allows ns to for¬ 
mulate several tests of NRQCD and the velocity scaling rules that it yields. 
In addition, we propose measurements of polarised gluon densities, Ag. The 
spin asymmetries have earlier been computed at lower orders for low trans¬ 
verse momentum processes and at high transverse momenta 0. 

The NRQCD factorisation formula for the inclusive production of heavy 
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quarkonium resonances H with 4-nionientuni P can be written as 



( 1 . 2 ) 


where $ is a flux factor. The coefficient functions Cij and Cij are computable 
in perturbative QCD and hence have an expansion in the strong coupling as 
(evaluated at the NRQCD cutoff). Although the matrix element is non- 
perturbative, it has a fixed scaling dimension in the quark velocity v. 

The fermion bilinear operators /Cj are built out of heavy quark fields 
sandwiching colour and spin matrices and the covariant derivative D. The 
composite labels i and j include the colour index a, the spin quantum number 
S, the number of derivative operators N, the orbital angular momentum L 
(found by coupling the derivatives), the total angular momentum J and the 
helicity J^. When the final state spin is unobserved, the hadron projection 
operator 


n(i?) = 


( 1 . 3 ) 


s 


(where s denotes hadron states with energy less than the NRQCD cutoff), is 
diagonal in these quantum numbers. Then it is clear that the operators /Cj 
and )Cj in eq. (B) are restricted to have equal L, S, J and J^, for both da 
and dAa. 

The J^-dependence of these matrix elements can be factored out using 
the Wigner-Eckart theorem— 


{K.n(/f)4) = 


( 1 . 4 ) 



where the factors of 1/(2J + 1) come from a Clebsch-Gordan coefficient. 
This is conventionally included in the coefficient function. We have used the 


notation of 0 for both diagonal and off-diagonal operators. The NRQCD 
power counting rule for these matrix elements is— 


d — 3 -|- -|- N' 2(E'(j -|- 


( 1 . 5 ) 
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where Ed and Md are the number of colour electric and magnetic transitions 
required to connect the hadronic state to the state /Ci|0). 

In section ^ of this paper we present the main results for Act after a 
telegraphic review of the threshold expansion technique m , the kinematics 
and the appropriate Taylor series expansion of the perturbative matrix ele¬ 
ment. This material is given in greater detail in our earlier paper [^. The 
phenomenology is discussed in the hnal section 


2 The Threshold Expansion 

We choose to construct the coefficient functions using the “threshold expan¬ 
sion” technique of This consists of calculating, in perturbative QCD, 
the matrix element Ai connecting the initial states to hnal states with a 
heavy quark-antiquark pair (QQ), and Taylor expanding the result in the 
relative momentum of the pair, q, after performing a non-relativistic reduc¬ 
tion of the Dirac spinors. The resulting expression is squared and matched to 
the NRQCD formula of eq. (|1.2|) by inserting a perturbative projector onto 
a non-relativistic QQ state between the two spinor bilinears. The coefficient 
of this matrix element is the required coefficient function. 

Symbolically— 


iEMhh'\M{h,h')\^ = E„ayK.n((30)q), 


( 2 . 1 ) 


where the left hand sides are Taylor expanded in q. Each factor of q Fourier 
transforms into a factor of the covariant derivative D on the right hand side. 
Since each matrix element on the right of eq. ( p.l|) corresponds to an unique 


matrix element in eq. (|1.2|) , the order up to which the Taylor expansion is 
to be performed is determined by the scaling of the non-perturbative matrix 
elements with v. Since we require a classihcation of operators by the angular 
momentum, it turns out to be very convenient to use the spherical tensor 
methods detailed in [0. 

In this paper we evaluate the polarised cross sections to order The 

Taylor expansion order, A -|- A' < 6 is obtained by setting d = 9 and Ed = 
Md = 0 in eq. (|1.5|). Furthermore, since we examine the leading term in 
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Figure 1: The ratios CqqjCgg (full line) and ACqq/Cgg (broken line) computed 
using the GRSV LO parton density set as a function of the CM energy 
y/S at zero rapidity. 


perturbation theory, the perturbative projector has only one term— 

n(QQ) = \QQ){QQ\. (2.2) 

In agreement with 0,0 we use the relativistic normalisation of states 

{Q{P,OQ{(1^V)\Q{P',OQ{(1',V')) = 4:EpEq{27r)^6^{p-p')S^{q-q'), (2.3) 

with the spinor normalisations = iq^rj = 1. Expanding Ep = Eq = 
y/m? + in allows us to write the spinor bilinears in terms of transition 
operators built out of the heavy quark held. The kinematics is very simple 
to leading order in as- The momenta of the initial particles are pi and p 2 - 
We take pi to he in the positive ^-direction and p 2 to be oppositely directed. 
The net momentum P = pi + p 2 . 

The 4-momenta of Q and Q {p and p respectively) are written as 

p = + Ljq^ and p = ^P — LjqP (2.4) 

Note that p^ = p^ = where m is the mass of the heavy quark. The 
space-like vector q is always dehned in the rest frame of the pair, and 
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boosts it to any frame. We shall use Greek indices for Lorentz tensors and 
Latin indices for Euclidean 3-tensors. 

In the next two sections we will write down the parton level cross sections, 
Ad for qq ^ H and gg ^ H (where iL is a quarkonium state). The hadronic 
cross sections, Act, are obtained by multiplying these by appropriate parton 
luminosities— 

A£gg = Ag{xi)Ag{x2), (2.5) 

/ 

where Xi = {2m/\/S) exp{y) and X 2 = {2m/y/S) exp(—y), y/S is the centre 
of mass (CM) energy of the colliding protons, y is the rapidity at which 
the quarkonium is observed, and Aqf{x), Aqf{x) and Ag{x) are the quark, 
antiquark (of flavour /) and gluon polarised densities. For unpolarised cross 
sections, a, we must similarly dehne unpolarised luminosities Cqq and Cgg 
using the corresponding unpolarised parton densities. As shown in Figure 
g, for > 20 GeV, Cgg « Cgg. Also, present day data indicates that 
•C Cgg. Consequently, the qq channel may be neglected for double 
polarised asymmetries to good precision. 


2.1 qq QQ 

We write down the polarised matrix element for the subprocess qq QQ for 
completeness, and because it provides a simple illustration of the techniques 
used. It is given exactly by the expression 


M 



M2 


X 


[v{P2, h')jgT^u{pi, h)] L>/ 

r . 4 . 1 

MQa^T°‘g - Q^Q{q ■ o')T°'g , 

L ' M + 2m^ ^ ’ 'J ’ 


( 2 . 6 ) 


where T“ is a colour generator, u and v are the light quark spinors with 
helicities h and h' respectively, and ^ and q are the heavy quark Pauli spinors. 
The equations of motion for the initial state quarks have been used to obtain 
the explicitly gauge invariant matrix element above. The desired Taylor series 
expansion is obtained by using the relation M^ = 4(m^ -|- g^) to expand all 
factors with M. 

The squared matrix element for this process is easily written down. The 
difference between the polarised and unpolarised cross sections follows when 
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one notices that the light quark helicities give an overall factor of 1 — hh' 
in the squared matrix element. As a result, Cij = —Cij and the parton 
level spin asymmetries are all —1. We list the parton level polarised cross 
sections— 




Ad-'" 

^ j/y 

qq 




9 

-5(s — 4m ) 


Sdm*^ 






^ 2 'pA’/’/Sr.O 3 q 2\ I lA^Ah/3p2^ 


54m‘^ 

1 f4 


S{s — 4m^) 


o|"ts“) + 






-vrcsi^s^,) 


1 jorcsf) + ^or(=B?) + 


where s is the parton CM energy. 



2.2 gg QQ 


The squared matrix element for the gg process is technically a little more 
complicated. We work in a class of ghost-free gauges called planar gauges 
[0- The density matrix for an initial state gluon of momentum p and helicity 
h in these gauges is given by 




~9^lv + 


1 . 

[pQ^u + PuVQj + 


p-V 


p-V 


-fil/pcr. 




( 2 . 8 ) 


The vector V dehnes the gauge choice. We write V = Cipi + C 2 P 2 , with 
C 1 /C 2 (9(1). We verify that all our results are gauge invariant by the 

explicit check that they do not depend on the arbitrary coefficients Ci and 
C2. 

The sum of the matrix elements arising from the three Feynman diagrams 
(s-channel gluon exchange, A4s, and t- and w-channel quark exchanges. Ait 
and Aiu) can be decomposed into three colour amplitudes— 


Ai — —g'^SabS + —g'^dabcD'^ + 7:9^ fabcF'^- 

b 2 2 


(2.9) 
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The colour amplitudes S and D involve only Ait + Aiu, whereas F involves 
Ms as well as Mt — Mu- 

In order to write down our results, we find it convenient to introduce the 
notation 

-4 = “ind S,j = A,Zj + AjZi - Bij + e, ■ tLzZiZj, 

( 2 . 10 ) 

where 


^ (^1 ■ ■ Pi - 62 ■ Li^i ■ P2) , 

Bij = 61 ■ Tje2 ■ Lj + 62 ■ TjCi ■ Lj. 


( 2 . 11 ) 


Here is the polarisation vector for the initial gluon of momentum pi. In 
order to identify all terms to order we need the colour amplitude S to 
order — 


5 = 


- (^) ^ i^^v) + (q^ecT^v) - (^) ^^mZn {q^q^^v) 


16 






M 




jmA-nA-p 


{q^q^q^i^a^p) 


- Az„,ZnZpZr {q^q'^q^q^^p) 


64 


'v ) 

^im. ^n. 


M 


^-im. Sn 


_j_ 2fYi^3m'^np 


ZrZs 


( 2 . 12 ) 

where ;§ is the unit vector in the z-direction. The amplitude D differs only 
through having colour octet matrix elements in place of the colour singlet 
ones shown above. For the colour amplitude F we need the expansion 


^ - (^) -4^™ (p"'i'T‘v) + (riXVV',,). (2.13) 

In all three colour amplitudes, the terms in A are spin singlet and those in 
S are spin triplet. 

The density matrix in eq. (|2.8|) yields 


III 

* 

FD 

FD 

<1 

i J 2 hh' hh'SjmSpjni 

= |Ho[^^?^]o - Ea=±2H![^^?^]-Ai 

III 

* 

1 C C* _ 

4 2 ^hh' Z>jm>~>j/yu' 

|[cTg][][cTtgt]0 + Yz\=±2[(^q]VF^q^]-x- 



(2.14) 
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We also find, A-A* = = 1/8 and ^-5* = A^-5* = 0. Conseqnently, 

the difference between the nnpolarised [^] and polarised cases lies solely in 
the flipped sign of the J = 2 part in eq. ( |2.14|) . Furthermore, even and odd 
terms in the three colour amplitudes of eq. (|2.9|) do not interfere with each 
other. 

The direct ^Si polarised subprocess cross section is 


A(rW(l) = ^[-eW(7)+|-ei''*(9) + -e^.''*(9) 


48 


16 


(2.15) 


where 


(p = 


4m2 


5(s — 4m^), 


(2.16) 


and denotes combinations of non-perturbative matrix elements from 

the colour amplitude a {= S, D or F) at order v'^ for the polarised subprocess 
cross section. These are given by 


= 


ei«(9) = 


2m? 

1 




5o°) + 






73) 


:V> 


2m^ . 
m* 


yi5( 


m" L 


(2.17) 


ei'*(9) = 


1 


^ (3 p2\ I ^ li) (3 

2m« is * ' ■' + 5 “ 


ph 


As a consequence of eq. (|2.14|) the coefficient of J = 2 matrix elements 
changes sign between the polarised and nnpolarised cases. The cross sections 
for ip', T and all other ^Si states need only the replacement of the appropriate 
matrix elements in eq. ( 2.171) . 

The ^Po polarised subprocess cross section is 


^^ggis) = ^ 


u^e?(5) + ie?(7) 
1 


+ (];|e|“(9) + (9) + (^e?(9) 


(2.18) 


























where combinations of non-perturbative matrix elements appearing in eq. 
(|2.18|) are 


ef{5) = 
ef{7) = 
ef{9) = 

ei%9) = 
ef{9) = 


2m‘^ 




4^3- 


nr 




1 


8m® 

1 


(3p3) + (3Po^®p5) 


9 


on^si) + 


2w? “ ' 2m4 L 


10^3 


3ore-Po') - jOfiri) 


5m‘^ 


orep^), 


18m® . 


0^°CSl) + 50fCD'i) 


3n2^ 


(2.19) 

The difference between these and the unpolarised cross sections of are only 
in the sign of the coefficient of the J = 2 components in the combinations 
0§°(9) and 0g’(9). 

For the ^Pi state we hnd 


^^ggis) = ^ 


Te«(9) + lex-(9) + lex-(9) 


( 2 . 20 ) 


where the combinations of non-perturbative matrix elements are, 


0F(9) = 
0^^(9) = 
Qf{9) = 


2m^ 

1 




Ofi^Sl) + 


1 


2w? “ ' 2m4 

ho?('q)+ ^ 


6m^ 


3m® 




( 2 . 21 ) 


ore si) 


+\oreDi)+\oreDi) 


The ^Po cross section is 


A(T»(s) = v? 


i4ey(5) + Te»(7) 

1 


+ 'l^er(9) + 4®”(8) + iiP(8) 


( 2 . 22 ) 
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where the combinations of non-perturbative matrix elements are 


0g^(5) = 

ef{7) = 

ef{9) = 
= 

OF (9) = 


5m^ 

2 


VTSm' 


orepi), 

:vrepi^pi), 


:orepi) - 


Zm* “ 75m» 

^ :Of (‘P/) + ^ 


(=>P|) + 


9 


2x/7 




3m® 


son^Po^) - -on^pi) 

5 

+ loreoi 


+\oreDi) - ?or(“p>i) 


(2.23) 

In the combination 0^(9), the coefficient of Clp(^P|) is different from the 
unpolarised in magnitude as well as sign, because it arises from several 
sources. Only the part from the combination AiS ■ S* flips sign. The other 
J = 2 terms change sign as expected. The J = 3 term in the F colour 
amplitude is also derived from the J = 2 part of the AS ■ S* term, and hence 
flips sign. 

The production cross section for ^5*0 states is— 


= 


7’ 


^e|(3)T ^0^.(5) 
1 


0^(7) + —0^(7) + —Ql(7) 
18 ^ 48 ^ 16 ^ 

where the combinations of non-perturbative matrix elements are 

1 


(2.24) 


0"s(3) = 
0^(5) = 
0^5(7) = 
0 "d(7) = 
0^(7) = 


2m^ 

1 




o'li's;;), 


PlCsi'si). 


3m® 

1 


4 

75' 


20U^si) + ^vn^sys^,) 


(2.25) 


oTCp^). 


6m^ 
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There is no J = 2 contribution to t]c to 0{v'^) and hence the parton level 
asymmetry is identically equal to 1 . 

The cross section for the ^Pi quarkonium state is— 






—0^(5) + —0^(7) 
48 ^ 48 ^ 

1 


+ 'l^§s( 9 ) + ^et(9)+4®F(9) 


(2.26) 


where the combinations of non-perturbative matrix elements is given by 


eys) 

ey 7 ) 

e|(9) 

ei,(9) 

§f(9) 


1 


2nP 

1 


VS' 


o5(‘sS). 


-,vicss, 'SI), 




2m? 

1 

3m® 


of('ss), 

20ICSI) + ‘S„‘) 


2m‘^ 


3oiepi) - joicpi' 


+ 


15m® 


O^C-Dl), 


6 m^ 




(2.27) 


The coefficient of the matrix element Og(^F 2 ) in Q%{9) flips sign as expected. 
However, the other J = 2 term, Og{^D 2 ), is unaffected as it comes from 


■ A*. 


2.3 jg QQ and 77 ^ QQ 

The matrix elements for the two processes yp —>■ QQ and 77 —>■ QQ are 
closely related to the gg amplitudes. It is easy to check that 




and 


= 




(2.28) 


where D and S are the colour amplitudes given in eq. (|2.9|) , and e is the 
charge of the heavy quark. 

The gg cross sections for the production of any quarkonium state can be 
obtained from those for the gg process, eqs. (|2.15| )- (p.27|) , by the following 
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prescription— replace in (p (eq. |2.16| ) by aas, delete the ©5 and Qp terms, 

and replace the colour factor 5/48 for the terms in ©d by 2. 

These computations may be applied to double polarised asymmetries in 
almost elastic photo-production of quarkonium states. Then the hadron level 
cross sections are obtained from these parton level quantities by multiplying 
the latter by the unpolarised or polarised gluon densities, g{x) and Ag{x) 
respectively. Here x = 2m ?is the proton mass and is the pho¬ 
ton energy in the rest frame of the proton. In phenomenological applications 
the diffractive parts of these cross sections have to be separated. Techniques 
have been developed at the HERA experiments to do this |]^ . 


The 77 cross sections are obtained with the prescription— replace aj 


m 


ip (eq. |2.1(j|) by delete the Qp and ©^ terms in eqs. ( p^.l5| )-( p.27| ), and 
replace the colour factor 1/18 for the terms in ©5 by 16. No parton densities 
enter for point-like photon cross sections. If, on the other hand, the cross 
sections are dominated by resolved photon processes, then the calculation is 
the same as for hadro-production. 


3 Phenomenology 

In spite of the large number of non-perturbative matrix elements in the hnal 
results, it is possible to make several quantitative and qualitative comments 
about the polarisation asymmetries. The most reliable predictions are ob¬ 
tained when the parton level asymmetries are ±1 (the ^S'o and ^Pi quarko¬ 
nium states) and the matrix elements drop out of the asymmetries. In these 
cases the polarisation asymmetries may be used to measure the polarised 
gluon densities. In several other cases large simplihcations occur. Some of 
the asymmetries can be used to test NRQCD or the scaling laws for matrix 
elements. 

Many of the matrix elements involved in the phenomenology are unknown. 
We develop a scaling argument which allows us to make rough estimates. A 
dimensional argument, neglecting possible logarithms of m and v, can be 
used to write 

(/C,n(i/)/C]) = RhY^^A^-v^, (3.1) 

where Dij is the mass dimension of the operator, d is the velocity scaling 
exponent in NRQCD (eq. El). A is the cutoff scale below which NRQCD 
is dehned (A ~ m), and Yij and Rh are dimensionless numbers. We have 
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Matrix Elements 

Values 

Rh 

08 '*(’Si) 

oR ('So) + ^oi'* (Sp„) 

(1.12 ±0.14) X 10-2 GeV^ 

(3.90 ± 1.14) X 10-2 GeV^ 

0.22 ±0.03 

0.17 ±0.05 

of es,) 

of ('So) + j^of (»Fo) 

(0.46 ± 0.08) X 10-2 GeV=^ 

(1.60 ±0.51) X 10-2 GeV^ 

0.09 ±0.02 

0.07 ±0.02 


Table 1: The values of various non-perturbative matrix elements in the non- 
relativistic normalisation El. These can be converted to the relativistic 
normalisation by multiplying by 4m. The dimensionless number Rh (eq. 
|Tl| ) has been extracted using A = m = 1.5 GeV, = 0.3, and F = 1 , 
except in the 2nd and 4th lines of the table where we used Y = 4.5. 


separated these two factors in order to specify that Yij contains all trivial 
numerical factors as well as Clebsch-Gordan coefficients, for example those 
which come from the approximate heavy-quark spin symmetry. Rh contains 
the irreducible minimum non-perturbative information. In the context of 
potential models, Rh would be related to some integral of the radial wave 
function of the state H. The scaling rules of NRQGD are a good guide to the 
importance of a given term if Rh depends mainly on H and is approximately 
independent of the operators /Cjj. 

In Table ^ we collect the values of some of the non-perturbative matrix 
elements quoted in the literature and extract from them the values of Rh- 
In the colour octet matrix elements were fitted to collider data. We 
have quoted the values obtained using the GRV LO parton densities. The 
good agreement between the values of Rh extracted from the two sets of 
octet matrix elements, for both H = J/ip and provide some basis for the 
working hypothesis that Rh is almost independent of the operator. In this 
section we shall accept it, keeping in mind that further checks are necessary. 

Assuming the constancy of Rh and using heavy-quark symmetry, we 
find that the non-perturbative matrix elements contribute approximately 
12R^w?v'^ to the Xi cross section and about (1 - 1 - -|- V2v^)R^rn?v^ to the 
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X 2 cross sections. Then we expect 


^(Xi) ^ 

cr(x2) l + v^ + 12v'^ 


0.45, 


(3.2) 


independent of y/S. This estimate^ is in reasonable agreement with the 
measured values in proton-nucleon collisions— 0.34 ± 0.16 at \/S = 38.8 
GeV [0 and 0.24 ± 0.28 at y/S = 19.4 []^. The measurements are also 
compatible with a lack of y/S dependence. In NRQCD this ratio cannot 
depend on the beam hadron. It turns out that the estimate in eq. (|3.2|) is 
not very far from the recently measured value in pion-nucleon collisions— 
0.57 ± 0.19 at y/S = 31.1 GeV |^. However, the experimental situation 


certainly needs clarihcation. 

This paper is concerned with low transverse momentum production. Our 
results are best suited to help design hxed target polarised experiments, for 
example by extracting the RHIG or HERA polarised beams for use on po¬ 
larised targets. The reason is that the polarisation asymmetries are likely 
to be largest at small y/S (Figure |[), yielding the best signals. Also, at 
smaller values of y/S the qq channel seems to be unimportant for the phe¬ 
nomenology, and the polarised gluon asymmetries can be studied directly. 
Since many aspects of low-energy quarkonium production can be tested by 
such experiments, it is worth considering them seriously. 

For quarkonium production, the parton level asymmetry in the gg 
channel is identically 1. For the yet to be discovered ^Pi resonances, this is 
true with corrections of order (from the term in eq. p.27|) . Only the 

colour amplitude S enters in 77 reactions; hence the polarisation asymmetry 
is exactly 1 for both these states. In yp collisions 


A — 
- 


Ap(x) 

9{x) 


[^Sq, ^Pi quarkonia). 


(3.3) 


For ^S'o states the above equation is exact, whereas there are order v'^ cor¬ 
rections for ^Pi states. 

^The number of terms allowed at higher orders in v can grow at most by 12 for each 
order. If the coefficient functions continue to be of order unity, then the all orders result 
for the Xi/X 2 ratio can be bounded by 0.67 for Xc and by 0.12 for Xb- The ratio in photo¬ 
production is close to unity in both these cases. In ||] an estimate of 0{a^) effects was 
used to show that the Xi/x 2 ratio could be about 0.3. 
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Figure 2: The asymmetries (a) ACgg/Cgg (full line) (b) {f^Cgg+ACqq)/{Cgg+ 
Cqq) (dashed line) and (c) (A£gg — ACqq)/{Cgg + Cqq) (dotted line) for zero 
rapidity as a function of the CM energy \fS using the GRSV LO parton 
density set |]TT[. The largest asymmetries are obtained at the lowest values 


of \fS, where all the asymmetries are dominated by the gluon asymmetry. 


For pp collisions, neglecting the qq channel for reasons explained earlier, 
we hnd that 


A — 

■^PP - 


1 + 0(v^ 


■'99 


c 


CSq, ^Pi quarkonia), (3.4) 


99 


where the order corrections come from the qq channel for the resonance, 
and in addition from the order n® matrix elements for the ^Pi resonance. 

The case of ^Pj states is very interesting phenomenologically. For yp — 
Xj reactions, we need to retain only the contribution from the colour ampli¬ 
tude D. The cross sections then begin at order n®, and we hnd that 

Ag{x) 


A — 


np 


^Pj quarkonia). 


(3.5) 


g(x) ej/ig) 

Since both the non-perturbative matrix elements involved are unknown, it 
might seem that there is no predictive power to the model. However, using 


heavy quark symmetry, and the scaling form of eq. (|3.1|) it appears that 
©d^(9) = 0 - As a result, we obtain the prediction 

,Ag{x) 


A^o = 


4x1 


— A^^ 

- 


[0 + O{v^)]- 


9{x) 


(3.6) 
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where the order terms come from the un-calculated order terms in the 
asymmetry and from corrections to heavy quark spin symmetry. 

We can make similar arguments for the polarisation asymmetries in 77 
collisions. We hnd that 


= l + 07 ) = -Af^. 


(3.7) 


The case of Xi is rnore complicated, since a ratio of two unknown matrix 
elements is involved. However, using heavy-quark spin symmetry this reduces 
to 


A^‘ = -- + 0(t? 


(3,8) 


77 y 

Unfortunately, no polarised 77 experiments have been planned to test this 
simple prediction of NRQCD and heavy quark spin symmetry. 

A straightforward application of the NRQCD scaling laws would lead us 
to the conclusion that the asymmetries for pp —> xop are given by 


4x0 

^PP 


\p 


AC 


gg 


C 


+ 0(v^ 


(3,9) 


gg 


where we have neglected the contribution of the qq channel. The asymmetry 
for xi production is 


= egH9) + eg(9) + eF(9) 
pp Qf{9) + ei^{9) + ef{9) 


AC 


gg 


C 


gg 


(3.10) 


The ratio of matrix elements can be estimated using heavy quark spin sym¬ 


metry and the scaling relations in eq. (|3.1| ). 0^'^(9) vanishes in this approx¬ 
imation and the terms 05 ‘^^( 9 ) come with opposite signs. The numerator is 
positive but small and we expect— 


4 x 1 


0.2 


C 


(3.11) 


gg 


The qq channel remains negligible even at '/S = 500 GeV. 

The ifg' asymmetry is straightforward. If the qq channel is neglected then 


4h' _ 

Wp - 


§0(7) 


©d(7) 


' D ' 

h'/ 


+ 0{v^), 


A^' = 

^IP 


eg (7) 

77 


^9{x] 

9{x) 0g(7) 


+ 0{v^), (3.12) 
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and the form of the non-perturbative matrix elements is the same as that 
for J/'^. The nnmerator can be shown to vanish using heavy quark spin 
symmetry and eq. (|3.1|) . A numerical estimate of the subleading order n® 
terms shows that the ratio of non-perturbative matrix elements is small (of 
the order of 0.05). At low '/S the asymmetry is dominated by the asymme¬ 
try in the gluon densities. However, at \/S of a few hundred GeV, the gluon 
asymmetry reduces to about 0.001, and the neglected qq channel becomes 
important. These terms appear with the opposite sign and reduce the asym¬ 
metry even further. As a consequence, we expect the xjj' asymmetry to be 
small for all \/S — 

0, Al;' ^ 0. (3.13) 


4^' 

^pp 


Any attempt to see structure in this asymmetry would probably require 
unreasonably large statistics. 

The J/ifj asymmetry seems to be enormously complicated because of the 
feed-down from radiative decays of the y states. However, a major sim- 
plihcation occurs because of the near vanishing asymmetry in direct J/'ip 
production (the argument is the same as for 'ip'). Thus the asymmetry comes 
entirely from the 20-40% of the cross section due to y decays. Taking into 
account the ratios of the production cross sections of y and the branching 
fractions for their decays into J/ip, we hnd that the xi ci-nd X 2 states con¬ 
tribute equally to J/pJ. Hence the J/pj polarisation asymmetry is expected 
to be approximately 


AJ/i’ 

pp 


-(0.15 ±0.05)-'^^^'' 


C 


(3.14) 


99 


We summarize the predictions made on the basis of the NRQCD scaling in 
eq. (|3.1|) and the assumption of Rh depending only on the hadron H — 


A^' -C 

^ ^pp ^ 


.AJ/i’ 

pp 


pp 


< A^° = 

PP 


AX2 — 

PP C„„ 1 


A^' 

7P 


4 Ah 
7P 


Axi 

7P 


^XO 

7P 


^X 2 

7P 


0 . 


(3.15) 


The asymmetry for J/ip in almost-elastic photo-production vanishes because 
the direct as well as feed-down contributions vanish. 

Alternative points of view can be taken about the scaling of the non- 
perturbative matrix elements. These proceed on the assumption that the 
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mass scale, A, appearing in the soft matrix elements are smaller than m. 
If the scale A in eq. ( 0 ) is taken to be the QCD scale Aqcd 0, then 
the relation v ~ as{mv) 0 implies that z for all m. 

The relative importance of terms then changes dramatically, since the low¬ 
est dimensional operators dominate. Since the dimension 4 terms 
appear in the y cross sections at order n® > zv^, they dominate the cross 
section. Heavy quark spin symmetry then gives the Xi/X 2 ratio to be 5/3, in 
agreement with the pion-nucleon observations. The same Xi/x 2 ratio is then 
predicted for charm and bottom mesons in hadro- and photo-production, 
and the result remains stable even when higher order terms are introduced. 
The rearrangement of the series also leads to ~ 7 ^ 0 and 

0^-^(9) ~ 0^‘^(9) 7 ^ 0. The consequences for the asymmetries are— 


Aii'7 





pp 

pp 

pp 

PP 

PP 

IP 


^XO 

IP 

~ AXi 

~ AX2 

~ ^7P 


AC 


99 


C 

^99 

Agjx) 

9{x) 


(3.16) 


and totally different from the NRQCD results summarized in eqs. am)- 
These differences make low-energy double polarised experiments a good test¬ 
ing ground for NRQCD, and the entire phenomenology of low-energy quarko- 
nium production. Other choices of the mass scale A have also been suggested 
[|T^ . These give somewhat different results. The asymmetries for ? 7 c,b and 

are independent of such assumptions, and hence form good measures of 
polarised gluon densities. 

In conclusion, low energy double polarised asymmetries are a good test¬ 
bed for understanding the origin of all observed systematics in fixed target 
hadro- and photo-production of charmonium. The high order computations 
presented here provide us with a set of processes which can be used to test 
aspects of NRQCD factorisation and scaling. In addition, and ^Pi char¬ 
monium and bottomonium states are well suited for measurements of po¬ 
larised gluon densities. Assuming standard NRQCD power counting laws, 
higher order terms are likely to be phenomenologically important for double 
polarisation asymmetries of ^Pi and excited ^Sq states in both hadro- and 
photo-production. 


18 







References 


[1] W.E. Caswell and G.P. Lepage, Phys. Lett., B 167 (1986) 437. 

[2] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev., D 51 (1995) 
1125; erratnm ibid., D 55 (1997) 5853. 

[3] E. Braaten, M. A. Doncheski, S. Fleming and M. Mangano, Phys. Lett., 
B 333 (1994) 548; 

D. P. Roy and K. Sridhar, Phys. Lett., B 339 (1994) 141; 

M. Gacciari and M. Greco, Phys. Rev. Lett., 73 (1994) 1586. 


[4] 

[5] 

[ 6 ] 
[7] 


S. Gnpta and K. Sridhar, Phys. Rev., D 54 (1996) 5545. 

M. Beneke and 1. Rothstein, Phys. Rev., D 54 (1996) 2005; 

M. Beneke, preprint GERN-TH/97-55, Piep-ph/9703429| . 

S. Gnpta and K. Sridhar, Phys. Rev., D 55 (1997) 2650. 

S. Gnpta and P. Mathews, preprint [hep-ph/9703370|, to appear in Phys. 


Rev., D. 

[8] S. Gnpta and P. Mathews, Phys. Rev., D 55 (1997) 7144. 

[9] O. Terayev and A. Tkabladze, preprint |hep-ph/ 9612301 . 

[10] E. Braaten and Y. Ghen, Phys. Rev., D 54 (1996) 3216. 

[11] M. Gliick, E. Reya and M. Stratmann and W. Vogelsang, Phys. Rev., 
D 53 (1996) 4775. 

[12] Yu. L. Dokshitzer, D. 1. Dyakonov and S. 1. Troyan, Phys. Rep., 58 
(1980) 269. 

[13] M. Beneke and M. Kramer, Phys. Rev., D 55 (1997) 5269. 

[14] K. Hagan-Ingram (E771 Gollaboration), unpublished; data quoted in 
the talk by M. Beneke p. 

[15] D. A. Bauer et ah, Phys. Rev. Lett., 54 (1985) 753. 

[16] V. Koreshev et ah, Phys. Rev. Lett., 77 (1996) 4294. 


19 











[17] M. Klein, preprint DESY 96-218, |hep-ph/9610398| ; 

B. Clerbaux, preprint DESY-96-162F; 

S. Schiek et al., HI Collaboration, preprint DESY-96-162E; 
S. Aid et al., HI Collaboration, Nucl. Phys., B 472 (1996) 3. 

[18] G. Schuler, preprint CERN-TH-97-012, |hep-ph/97U223(J| . 


20 






